Emergence of the N=16 shell gap in 21 
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The spectroscopy of 21 has been investigated using a radioactive 20 O beam and the (d,p) reaction 
in inverse kinematics. The ground and first excited states have been determined to be J 7r =5/2 + and 
J n = l/2 + respectively. Two neutron unbound states were observed at excitation energies of 4.76 ± 
0.10 and 6.16 ± 0.11. The spectroscopic factor deduced for the lower of these interpreted as a 3/2 + 
level, reveals a rather pure 0d 3 / 2 single-particle configuration. The large energy difference between 
the 3/2 + and l/2 + states is indicative of the emergence of the N=16 magic number. For the higher 
lying resonance, which has a character consistent with a spin-parity assignment of 3/2 + or 7/2 _ , a 
71% branching ratio to the first 2 + state in 20 O has been observed. The results are compared with 
new shell model calculations. 

PACS numbers: 21.10.Hw 21.10.Jx 23.20.Lv 25.60.Je 27.30. +t 29.38. Gj 



The magic numbers that explain the structure of nu- 
clei close to stability have their origin in the gaps created 
by the single-particle eigenstates of the mean-field. The 
structure of exotic nuclei is now known to often differ 
from that near stability, exhibiting an evolution of the 
shell closures or even a reordering of the levels [l|. The 
single-particle properties can be substantially modified in 
light neutron-rich nuclei by, in part, to the combined ac- 
tion of the central component and the tensor part of the 
effective nucleon-nucleon (NN) interaction 0, Fur- 
thermore, near the drip line, other effects become im- 
portant including many-body correlations [irQ], an d the 
influence of the scattering to the continuum in weakly 
bound and unbound states 0-@|- 

The neutron-rich oxygen isotopes represent an intrigu- 
ing means to study the interplay of such effects. The last 
bound oxygen isotope is 24 0, which reinforces N=16 as 
a shell gap, but the addition of a single proton moves 
the neutron-drip line for the fluorine isotopes to at least 
31 F (N=22) [lO] . The N=16 gap is produced by an in- 
crease in the spacing between the lsi/ 2 and the 0d 3 / 2 
neutron orbitals. A precise knowledge of how these or- 
bitals evolve is crucial in order to predict the position 
of the neutron-drip line. Although fundamental in de- 
termining the binding energy of 26 < 28 0, the energy of the 



v0d 3 /2 single-particle orbital is still the object of disagree- 
ment between different shell model interactions. Several 
phenomenological interactions, such as the USDB [l2[ 
and SDPF-M are able to successfully predict 24 as 
the last bound oxygen isotope. However, when derived 
from microscopic nucleon-nucleon (NN) forces, the posi- 
tion of the drip line is shifted to 28 O [4j. The inclusion 
of three-nucleon forces (3NFs) is known to influence the 
location of the drip line as pointed out first in ab-initio 
calculations for light nuclei by Hagen et al. [5J and, more 
recently by Otsuka and collaborators The influence 
of continuum coupling, while being important in describ- 
ing the low-lying states in 24 0, does not appear to affect 
the position of the drip line [§]. 

Previous work investigating the i/0c? 3 / 2 orbital em- 
ployed transfer of a neutron onto 22 O [14|. The first 
3/2+ level in 23 was found at 4.00(2) MeV consistent 
with shell model calculations using the SDPF-M [l3| and 
USDA [12] interactions. In 25 0, the location of the un- 
bound v0djU2 ground state, observed in proton knockout 
from 26 F 11511 could only be reproduced by the USD inter- 
action [ll| which in contrast does not predict correctly 
the neutron-drip line. In the present work we have em- 
ployed the d( 20 O,p) reaction as a means to shed further 
light on the evolution of the vQd 3 / 2 strength. In the naive 



2 



shell model picture 20 O has 4 neutrons in the vOdh/2 or- 
bital. The stripping of a nucleon from the target is thus 
expected to populate the v0d 5 / 2l v\s\ii and ^0d 3 / 2 or- 
bitals. The size of the N=16 shell gap is related to the 
energy difference between the 3/2 + and l/2 + states. Ex- 



isting data on O [TJj, [18| provide evidence for a low- 
lying 3/2+ state, but this is expected to include signifi- 
cant core excitation rather than a strong ^0^3/2 single- 
particle structure. We report here on the measurement of 
the observation of a higher lying 3/2 + level that carries 
significant v0d 3 /2 strength. 
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FIG. 1: Upper: Barrel excitation energy spectrum of the un- 
bound states in 21 from 20 O-|-p coincidences. Protons were 
detected between [95-145] degrees. Bottom: Events corre- 
sponding to protons in coincidence with 1684 keV deexcita- 
tion 7 rays from 20 O. The inset of the bottom figure shows 



the 7-ray spectrum gated on 
the TIARA Si array. 
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The experiment was performed at the SPIRAL facility 
of GANIL in Caen using a 10.53 AMeV 20 O beam with 
an average intensity of 10 4 pps. The beam impinged on a 
0.59 mg/cm 2 thick CD 2 target. A detailed description of 
the experimental set-up can be found in references 
2 If . The (d,p) reaction channel was selected using the 
identification of the beam-like residue using the VAMOS 
spectrometer and the kinematics of the proton in the 
position sensitive silicon array TIARA [19j which sur- 
rounded (37°-145°) the target. The excitation energy 
in 21 and the centre-of-mass (cm) scattering angle was 
obtained from the measured energy and laboratory angle 
of the recoiling protons. The bound first-excited state of 
21 was measured at 1.21 ± 0.04 MeV. As in our earlier 
work [IH21I, the target was also surrounded at 90° by 
4 EXOGAM clover detectors in a close packed geometry 



(efficiency was 10% at 1 MeV). The coincident gamma- 
energy spectrum provided a measurement of the excita- 
tion energy of 1.213 ± 0.007 MeV. The values obtained 
both via the detection of the protons and the gamma ray 
compare well with previously established energy of 1.218 
± 0.004 [l8j]. Any other bound states populated were 
over an order of magnitude weaker. 

The excitation energy spectrum deduced for states ly- 
ing above the neutron emission threshold (Sn=3.806 ± 
0.012 MeV) in 21 is shown in the upper part of Fig. 
[U The solid line represents the sum of two resonances 
(dashed and pointed lines) and a contribution from di- 
rect break-up ( 20 O+n+p). The resonance line-shape em- 
ployed was a Breit-Wigner convolved with a Gaussian 
of fixed width (FWHM exp ) that represented the experi- 
mental resolution accounting for the beam characteristics 
(spot size and energy spread), the energy and angular 
straggling in the target, detector resolutions and kine- 
matical effects. The three-body contribution, taken here 
to mean a background from the breakup of the deuteron 
( 20 O+n+p) , was obtained by uniformly sampling the 
phase-space of such a decay. The normalisation of this 
contribution was a free parameter in the fit. A Monte- 
Carlo simulation [l9| based on Geant4 was performed to 
determine the experimental resolution. The simulations 
were validated by comparison with the measured reso- 
lution of the state at 1.21 MeV, FWHM e:E p=0.77 ±0.17 
MeV and FWHM sim =0.75 MeV. The simulations were 
then used to estimate the resolution at higher excitation 
energies. The adjustement to the data yielded E reSl = 
0.960 ± 0.10 MeV, Ti=0.560 ± 0.20, and E reS2 =2M 
± 0.11 r 2 =0.390 ± 0.44. The corresponding excitation 
energies, as listed in Table HI of the unbound states are 
4.76 ± 0.10 and 6.16 ± 0.11 MeV. The bottom part of 
Figfflshows the excitation energy spectrum for the events 
for 7-particle-fragment triple coincidences, ( 20 O+p+7). 
Taking into account the 7-ray efficiency, the second res- 
onance was found to decay to the 1.684 MeV (2 + ) first 
excited state in 20 O with a branching ratio, r(2 + )/r tot , 
of 0.71 ±0.22. 

Shell model calculations presented in table El sug- 
gest that the states expected to be populated above 
the neutron threshold include the 3/2 + and, possibly, 
negative parity states from the /p-shell. Any £ = 1 
single-particle resonances at the measured excitation en- 
ergies are expected to be very broad (r sp (4.76)=3.4 MeV, 
r sp (5.5)=100 MeV). Therefore, the measured widths of 
the observed resonances already excludes them from be- 
ing £ = 1 states with a reasonable spectroscopic factor. 

As in reference [2(| , elastic scattering was used to cal- 
ibrate the beam current and number of deuterons in 
the target. The optical- model parameters obtained from 
d± 26 Mg at 12 AMeV ^ were used and the final uncer- 
tainty in the normalisation was estimtated to be at most 
2%. 

Figure [2] displays the angular distributions for the pro- 
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FIG. 2: Differential cross sections for the bound a) and b) 
and unbound c) and d) states in 21 compared to ADWA 
calculations. For the adopted f-transfers see the text and 
Table I. Only statistical uncertainties are included. 



TABLE I: Results from the present work for the states ob- 
served in 21 0. BR=r(2+)/T tot for 20 O+n decay. 



E x (keV) 


t 


J" 


BR 


C 2 S 





2 


5/2+ 




0.34 ± 0.03 


1213 ± 7 





1/2+ 




0.77 ± 0.09 


4760 ± 100 


2 


3/2+ 




0.58 ± 0.06 


6160 ± 110 


2 


3/2+ 


0.71±0.22 


0.30 ± 0.05 




3 


7/2" 


0.71±0.22 


0.20 ± 0.02 



tons leading to the bound and unbound states in 21 
compared to reaction model calculations. The theoret- 
ical cross sections were computed within the Adiabatic 
Distorted Wave Approximation (ADWA) using an adi- 
abatic potential for the d+ 20 O entrance channel [26| . 
The Chapel-Hill parameterisation (CH89) [27[ was used 
to describe the nucleon-nucleus potentials for both the 
entrance and outgoing channels. All the calculations for 
the unbound states were performed using the Vincent 
and Fortune prescription [32| implemented in DWUCK4 
[33| . The neutron single-particle form factors were ob- 
tained using a Woods-Saxon potential with the depth 
adjusted to reproduce the experimental binding energy 
of each level. Checks on the results for the unbound 
states were performed using bins in the continuum repre- 
senting each resonance and implemented in the FRESCO 
code [34{ . Consistent results for the angular distributions 
were obtained with both approaches. The measured ex- 
citation energies, spectroscopic factors, and the assigned 
spins and parities are listed in Table I. 



The angular distribution for the transfer to the ground 
state is well reproduced by an I = 2 transfer, in agree- 
ment with a 5/2+ assignment made in earlier studies 
17 . llj]. The differential cross sections for the 1.21 MeV 
state are only consistent with an £ — transfer, which 
permits a definite l/2 + assignment to this state. The 
spectroscopic factors were obtained by minimising the 
X 2 for all the data points. The spectroscopic factors of 
the ground and first excited states were determined to 
be 0.34 ± 0.03 and 0.77 ± 0.09, respectively, in agree- 
ment with the predictions using the USDA, SDPF-M and 
WBP interactions (Table|Il]). The 5/2+ and 1/2+ bound 
states carry most of the available strength of the 0d$/2 
and \si/2 orbitals, respectively, based on the vacancies 
of the projectile. The occupancy of the lsi/2 orbital is 
smaller than unity as would be expected in a naive shell 
model picture. This is consistent with this orbital be- 
ing partially occupied in the projectile wave function as 
shown by a complementary study of the 20 O(d,t) 19 O re- 
action [301 ] . 

The angular distribution for the resonant state at 
Ej;=4.76 MeV is best reproduced by an I — 2 transi- 
tion, leading to a spin-parity assignment of </ 7r = (3/2+, 
5/2+). Comparisons with shell model calculations using 
different interactions predict no significant spectroscopic 
factor strength for a 5/2+ state in this energy region (Ta- 
ble |n|). States with negative parity /p-configurations are 
predicted to lie 1 MeV higher in excitation energy than 
the strong, ^0d 3/2 , 3/2+ state by the WBP and SDPF- 
M interactions. Therefore, it is concluded that the state 
at 4.76 MeV is very likely 3/2+. The corresponding ex- 
perimental spectroscopic factor of 0.58 ± 0.06 for I — 2 
is in very good agreement with the USDA value of 0.68, 
while, as listed in Table [Til the WBP interaction pre- 
dicts this state to carry a somewhat larger fraction of 
the 0c?3/2 single-particle strength. The SDPF-M interac- 
tion predicts the corresponding state to have a very low 
spectroscopic factor and with the strength concentrated 
in a higher lying state. 

For the second unbound state at 6.16 MeV, the an- 
gular distributions for I — 2 and 3 transfer agree with 
the data equally. The corresponding spectroscopic fac- 
tors are 0.30±0.05 and 0.20±0.02. For I = 1, no reso- 
nance can be generated this far above the threshold, as 
such, a very small spectroscopic factor would, therefore, 
be needed for the width of an I = 1 state to be consis- 
tent with the observed width of 0.39 MeV. As the decay 
of the 6.16 MeV state to the ground state represents only 
0.30 ± 0.22 of the decay, an £ = 1, (3/2") assignment is 
very unlikely. From Table ITfl the states expected to be 
strongly populated in the (d,p) reaction can be identified. 
Leaving aside the 3/2 2 + (3/2 3 + state in SDPF-M) iden- 
tified with the 4.76 MeV state, the remaining candidates 
are the other 3/2+ states together with the 7/2~ and 
3/2~ states. For all of these, the energy gap is 1.5 MeV 
relative to the 3/22+. On the other hand the only states 
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TABLE II: Shell Model calculations for 21 using the USDA/SDPF-M/WBP interactions for the sd-pf shell states. 



J" (nlj) 


i 


USDA [12] 
£4keV) C 2 ^ 


SDPF-M [13] 
^(keV) C 2 S 


WBP [31] 

^(keV) C 2 S 


l/2 + lSl/2 


1 

2 


1277 
6981 


0.83 
0.00 


1451 
6695 


0.82 
0.00 


1331 

5988 


0.81 
0.00 


3/2+ 0d 3 /2 


1 

2 
3 
4 


2186 
5522 
5671 
8195 


0.00 
0.68 
0.13 
0.11 


1926 
5446 
6008 
8666 


0.00 
0.08 
0.77 
0.09 


2189 
4836 
5796 
7897 


0.00 
0.82 
0.03 
0.06 


5/2+ 0d 5 /2 


1 

2 
3 
4 




3129 
4856 
7537 


0.34 
0.04 
0.00 
0.01 




3034 
4796 
7296 


0.33 
0.04 
0.00 
0.00 


0.0 

3149 
4965 
6547 


0.34 
0.04 
0.00 
0.01 


3/2" lp3 /2 


1 
2 
3 






7003 
8427 
8760 


0.48 
0.15 
0.00 


6223 
6905 
7481 


0.36 
0.10 
0.02 


7/2- l/ 7/2 


1 

2 






6812 
8698 


0.70 
0.11 


5942 
6469 


0.12 
0.05 



in Table [TTI favouring decay to the 20 O(2 + ) turn out to 
be the 3/24 + , the 7/22 _ and the 3/2^~ states. For the 
3/24 + , for example, the three interactions predict values 
of 0.69, 0.74 and 0.76 for r(2+)/r tot , but are associated 
with a spectroscopic factor of only 0.1 (Table HI|) . It 
seems, therefore, that the shell model does not mix the 
various configurations in the 3/2 + states correctly. As- 
suming the state at 6.16 MeV to be J 71 ' = 3/2 + , then 
the 0d 3 / 2 strength would be split between the 4.76 and 
6.16 MeV states in a manner that is not predicted by the 
calculations. The other possibility is a 7/2~ assigment, 
with a wavefunction carrying a small part of the 
strength and dominated by core excitations coupled to a 
neutron. 

As noted earlier, an estimate of the limit of the size 
of the N=16 shell gap in 21 can be deduced from the 
difference in excitation energies of the l/2 + (Ex=1.21 
MeV) and the 3/2+ (Ex=4.76 MeV) states. The value 
here of 3.5 MeV is close to that of 4 MeV in 23 O [l4|. 
The invariant mass spectra of the unbound states in 23 O 
populated in the d( 22 0,p) reaction also shows a second 
resonance at roughly 1.3 MeV distance to the first one 
compare to ours at 1.4 MeV. A 3/2~ assignment was 
proposed for this state [3]. However, the width of the 
resonance is too narrow to be in agreement with the spec- 
troscopic factor of 1 deduced from the magnitude of the 
cross section. This inconsistency might arise from the 
weakly bound approximation used for the final state in 
the reaction calculations. 

Summarising, we have presented results for neutron 
transfer populating the bound and unbound states in 
21 0. The first excited state was assigned spin parity l/2 + 
and the £ = 2 transfer to the ground state supports a 
5/2+ assignment. The spectroscopic factors derived for 
both states are in good agreement with the shell-model 
predictions. Two neutron unbound states have been ob- 



served at 4.76 and 6.16 MeV. The 4.76 MeV state is 
identified as 3/2 + and, as such, carries 58% of the 
v0d 3 / 2 strength. The state at 6.16 MeV corresponds ei- 
ther to a large fraction of the missing ^0d 3 / 2 strength 
or to an intruder /-state from the p/-shell. The shell 
model is unable to provide an adequate explanation for 
this level. This poses a challenge for the shell model 
and the interactions used to describe the evolution of nu- 
clear structure in this region. Ideally the issues discussed 
here should be investigated using less phenomenological 
models that include explicitely the physics of interest, 
such as the continuum. The triple-coincidence technique 
has been shown to be a very useful tool to probe the 
structure of unbound states with sizeable core excitations 
components. This suggests the equipment being devel- 
oped at emerging ISOL facilities for studying transfer re- 
actions should incorporate highly efficient 7-ray arrays, 
even when studying unbound states. 
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